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ABSTRACT 

Transient black hole candidates are interesting objects to study in X-rays as these sources show rapid evo¬ 
lutions in their spectral and temporal properties. In this paper, we study the spectral properties of the Galactic 
transient X-ray binary MAXI J1659-152 during its very first outburst after discovery with the archival data of 
RXTE Proportional Counter Array instruments. We make a detailed study of the evolution of accretion flow 
dynamics during its 2010 outburst through spectral analysis using the Chakrabarti-Titarchuk two-component 
advective flow (TCAF) model as an additive table model in XSPEC. Accretion flow parameters (Keplerian 
disk and sub-Keplerian halo rates, shock location and shock strength) are extracted from our spectral fits with 
TCAF. We studied variations of these fit parameters during the entire outburst as it passed through three spectral 
classes: hard, hard-intermediate, and soft-intermediate. We compared our TCAF fitted results with standard 
combined disk black body (DBB) and power-law (PL) model fitted results and found that variations of disk 
rate with DBB flux and halo rate with PL flux are generally similar in nature. There appears to be an absence 
of the soft state unlike what is seen in other similar sources. 

Subject headings: X-Rays:binaries - stars individual: (MAXI J1659-152) - stars:black holes - accretion, ac¬ 
cretion disks - shock waves - radiation:dynamics 


1. INTRODUCTION 

Compact objects, such as black holes (BHs) and neutron 
stars are identified by electromagnetic radiations emitted from 
the accreting disk formed by matter from their companion. 
This matter is accreted either through Roche-lobe overflow 
or winds. Some of these objects are transient stellar mass 
X-ray binaries with a low mass star acting as a donor star. 
Study of these objects in X-rays is very interesting as they un¬ 
dergo rapid evolution in their timing and spectral properties 
which are strongly correlated to each other. There are a large 
number of articles by several groups (see, e.g., McClintock & 
Remillard, 2006; Belloni et al. 2005; Nandi et al. 2012; Deb¬ 
nath et al. 2013, etc.) which discuss variation of spectral and 
temporal properties of these transient black hole candidates 
(BHCs) during their X-ray outbursts. In general, it has been 
found that these objects show different spectral states (hard, 
hard-intermediate, soft-intermediate, soft, etc.; McClintock & 
Remillard, 2006) and low and high frequency quasi-periodic 
oscillations (QPOs) in power-density spectra (PDS) in some 
of these spectral states (Remillard & McClintock, 2006). Dif¬ 
ferent branches of X-ray color-color and Hardness-Intensity 
Diagrams (HIDs; Maccarone & Coppi, 2003; Belloni et al. 
2005; Debnath et al. 2008, etc.) are related to different spec¬ 
tral states of the outburst phases. It has also been observed 
by several authors (see, Mandal & Chakrabarti, 2010; Nandi 
et al., 2012; Debnath et al., 2013 and references therein) that 
these observed spectral states show hysteresis loops during 
their spectral evolutions of an entire epoch of the outburst of 
these transient BHCs. Simple model fits of accretion rates 
using two component advective flow (TCAF) solution of 
Chakrabarti-Titarchuk (1995, hereafter CT95) by Mandal 
and Chakrabarti (2010) indicated that indeed the accretion 
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rates vary differently in the rising and declining states. 

In the literature, there are a large number of theoretical 
or phenomenological models, which describe accretion flow 
dynamics around a BH. It is well known that emitted ra¬ 
diation contains both thermal and non-thermal components. 
The thermal component is multi-color blackbody type emit¬ 
ted the standard Keplerian disk (Shakura & Sunyaev, 1973; 
Novikov & Thorne, 1973) and the other is a power-law com¬ 
ponent, originated from the so-called ‘Compton’ cloud (Sun¬ 
yaev & Titarchuk, 1980; 1985). This component is composed 
of hot electrons and is cooled down by repeated Compton 
scatterings of the low energy (soft) photons. There are many 
speculations about the nature of this Compton cloud which 
range from it being a magnetic corona (Galeev, Rosner & 
Viana, 1979) to a hot gas corona over the disk (Haardt & 
Maraschi, 1991; Zdziarski et al., 2003). CT95, in their TCAF 
solution, considers that the CENtrifugal pressure supported 
BOundary Layer or CENBOL plays the role of the Compton 
cloud. This CENBOL happens to be the post-shock region of 
the low-angular momentum halo in which a standard Keple¬ 
rian disk remains immersed while emitting soft photons. The 
shock in the halo forms due to the piling up of matter behind 
the centrifugal barrier of the low angular momentum accre¬ 
tion flow component having a sub-critical viscosity parameter 
(Chakrabarti 1990ab, hereafter C90ab, 1996). These shocks 
are found to be stable even under non-axisymmetric pertur¬ 
bations (Okuda, Teresi & Molteni, 2007). The other compo¬ 
nent of TCAF is an optically thick standard (SS73) Keplerian 
(disk) component which is formed in flows with super-critical 
viscosity parameter (C90ab). Of course, this disk also has 
to pass through the inner sonic point to satisfy the boundary 
conditions on the black hole horizon (C90ab; see also, Mu- 
chotrzeb & Paczynski, 1982). Formation of TCAF from a 
in single simulation (see, Giri & Chakrabarti, 2013; Giri, Garain 
& Chakrabarti, 2015, and references therein) show that it has 
a stable configuration. Also, in Mondal et al. (2014a) a self- 
consistent transonic solution of TCAF in presence of both 
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cooling and outflows is obtained. 

Recently, after the inclusion of this TCAF solution (CT95; 
Chakrabarti, 1997, hereafter C97) in HEASARC’s spec¬ 
tral analysis software package XSPEC as an additive ta¬ 
ble model Debnath, Chakrabarti & Mondal (2014, hereafter 
DCM14), Debnath, Mondal & Chakrabarti (2015, hereafter 
DMC15) and Mondal, Debnath & Chakrabarti (2014b, here¬ 
after MDC14) obtained a clearer picture about the accretion 
flow dynamics around BHCs as they find evidences of system¬ 
atically varying accretion rates of the standard disk and the 
halo and the shock location on a daily basis. From the TCAF 
fitted spectrum, one can obtain two accretion rates (disk and 
halo), shock locations and shock strength. These parameters 
also give us information about the frequency of QPOs. Tran¬ 
sitions of various spectral states which are observed during 
the outburst phases of a transient BHC can be identified by 
special behaviour of accretion rate ratio (ARR) and the nature 
of observed QPOs. 

Newly discovered MAXI J1659-152 is an interesting black 
hole binary to study because it is the shortest orbital period 
BHC observed till date (Kuulkers et al„ 2010; 2013). The 
source was first observed by MAXI/GSC instrument on 25th 
Sept. 2010 at the sky location of R.A. = 16 /! 59"' 10 s , Dec 
= -15°16'05" (Negoro et ah, 2010). The source was simulta¬ 
neously observed by SWIFT/BAT instrument roughly at 17° 
above the Galactic plane (Mangano et al. 2010). Kalamkar 
et ah (2011) defined the source as a BHC, based on their 
combined optical and X-ray spectral study, which was ini¬ 
tially thought to be a Gamma-Ray Burst and was named as 
GRB100925A). Kuulkers et ah (2013) based on their de¬ 
tailed study of the X-ray intensity variation of observed ab¬ 
sorption dips (Kennea et ah 2010), confirmed MAXI J1659- 
152 to be a short orbital period black hole binary of period 
= 2.414 + 0.005 hrs. 

MAXI J1659-152 showed X-ray flaring activity in 2010, 
other than low-level activity in 2011 which continued for 
~ 9 months. During this period, the source was extensively 
studied in multi-wave band, such as various X-ray (Munoz- 
Darias et ah 2011; Kalamkar et ah 2011; Yamaoka et ah 
2012; Kuulkers et ah 2013), optical/IR (Russel et ah 2010; 
Kaur et ah 2012), and radio observatories (Miller-Iones et 
ah 2011; Paragi et ah 2013). van der Horst et ah (2013) 
made multi-band campaign to explore multi-wavelength prop¬ 
erties of the source during this outburst. Physical parameters, 
such as distance, disk inclination angle, masses of the source 
and the companion are estimated using various methods. The 
most acceptable ranges of distance and disk inclination angle 
are 5.3 - 8.6 kpc and 60 - 80 degree (Yamaoka et ah 2012; 
Kuulkers et ah 2013) respectively. Although Shaposhnikov et 
ah (2011) predicted the mass of the BH ( Mbh ) as 20 + 3 M 0 , 
the preferable range of mass of the source and companion are 
3-8 M 0 (Yamaoka et ah 2012) and 0.15- 0.25 M e (Kuulkers 
et ah 2013) respectively. In this paper, we use mass of the BH 
as 6 M 0 . 

We study spectral and timing properties of the source dur¬ 
ing its 2010 main outburst phase, which continued for ~ 
1.5 months, using RXTE PCA archival data. Temporal prop¬ 
erties of the BHC along with the evolution of QPO frequency 
during declining phase of the outburst are presented in Molla 
et ah (2015, hereafter Paper II). 

The paper is organized in the following way: in the next 
Section, we briefly discuss observation and data analysis pro¬ 
cedures using HEASARC’s HeaSoft software package. In §3, 
we present results of spectral analysis using TCAF fits file 


as an additive table model in XSPEC and variation of differ¬ 
ent flow parameters extracted from model fits. Here, we also 
compare combined DBB and PL model fitted spectral analysis 
results with that of the TCAF fitted analysis results. Finally, 
in §4, we present a brief discussion and make our concluding 
remarks. 

2. OBSERVATION AND DATA ANALYSIS 

We analyze the data of 30 observational IDs starting from 
the first day of RXTE PCA observation, namely, 2010 
September 28 (Modified lulian Day, i.e., MJD = 55467) to 
2010 November 11 (MJD = 55508). Data reduction and anal¬ 
ysis are done using HEASARC’s software package HeaSoft 
version HEADAS 6.15 and XSPEC version 12.8. To analyze 
archival data of the RXTE PCA instrument, we follow the 
standard data analysis techniques as done by Debnath et ah 
(2013,2015). 

For spectral analysis, Standard2 mode Science Data of PCA 
(FS4a*.gz) are used. Spectra are extracted from all the layers 
of the PCU2 for 128 channels (without any binning/grouping 
the channels). We exclude HEXTE data from our analysis, as 
we find strong residuals (line features) in the HEXTE spec¬ 
tra at different energies. This could be due to the fact that 
the ‘rocking’ mechanism for HEXTE stopped. So, we re¬ 
strict our spectral analysis with the PCA data for the energy 
range of 2.5 - 25 keV only. In the entire PCA data analysis, 
we include the dead-time correction and also the PCA break¬ 
down correction (because of the leakage of propane layers of 
PCUs). The “runpcabackest” task was used to estimate the 
PCA background using the latest bright-source background 
model. We also incorporated the pcasaa-history file to take 
care of the SAA data. To generate the response files, we used 
the “pcarsp” task. Detailed analysis will be discussed in Paper 
II. 

The 2.5-25 keV PCA background subtracted spectra are fit¬ 
ted with TCAF based model fits file and with combined DBB 
and PL model components in XSPEC. Individual flux contri¬ 
butions for the DBB and PL model components are obtained 
by using the convolution model ‘cflux’ technique. For the en¬ 
tire outburst, we keep hydrogen column density (N#) fixed at 
3.0xl0 21 atoms cm ~ (Munoz-Darias et al., 2011) for absorp¬ 
tion model wabs. We also assume a fixed 1.0% systematic 
instrumental error for the spectral study during entire phase 
of the outburst. After achieving the best fit based on reduced 
chi-square value ix 2 red ~ 1), ‘err’ command is used to find 
90% confidence error values for the model fit parameters. In 
Appendix, Table I, we mention average values of these two 
errors in superscript. 

For a spectral fit, using the TCAF based model, one needs 
to supply five model input parameters, other than the normal¬ 
ization constant. These parameters are: i) black hole mass 
(Mbh) in solar mass (M 0 ) unit, if) sub-Keplerian rate (ni/, in 
M Edd), Hi) Keplerian rate (m d in Eddington rate M Edd), iv) lo¬ 
cation of the shock (X s in Schwarzschild radius r g =2GM/c 2 ), 
v) compression ratio (R) of the shock. The model normal- 
• R 2 ? 

ization value (norm) is jfi^sin(i), where ‘/C’ represents an 
effective area of the emitting region (purely on dimensional 
ground), D is the source distance in 10 kpc unit and i is the 
disk inclination angle. In order to fit a black hole spectrum 
with the TCAF model in XSPEC, we generate model fits file 
(TCAFO.l.fits) using theoretical spectra generating software 
by varying five input parameters in CT95 code. We then in¬ 
clude it in XSPEC as a local additive model. A brief dis- 
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cussion of the TCAF model, its present development and a 
detailed description of the range of input parameters and gen¬ 
eration procedure of the current version (vO.l) of the TCAF 
fits file are given in DCM14 and DMC15. For the spectral 
analysis with TCAF, mass of the black hole is frozen at 6 M 0 . 

3. RESULTS 

Accretion flow dynamics of a transient BHC can be well un¬ 
derstood by model analysis of spectral and temporal behaviors 
of the source during its outburst phase. Here, we present re¬ 
sults of spectral analysis based on TCAF and compare with 
combined DBB and PL model fitted results. A combined 
DBB and PL model fitted spectral analysis, though fitted well 
throughout and in some cases better than TCAF, only give 
gross properties of the disk such as fluxes from different com¬ 
ponents. However, TCAF goes one step further in extract¬ 
ing the detailed flow parameters, such as two disk rates and 
shock properties. Furthermore, transitions of spectral states 
are more conspicuous in terms of the fitted parameters. Thus, 
to study accretion dynamics around BHCs, there appears to 
be certain definite advantages in fitting with TCAF solution. 
A shortcoming of TCAF fit with the current version (vO.l) 
is that as the spectra become softer, the fit tends to worsen, 
mainly indicating that the importance of the halo component 
is reducing. In our next version this would be taken care of by 
self-consistently cooling down the second component in order 
that the flow automatically tends to have a single component. 

All 30 observational IDs spreaded over the entire period of 
the 2010 outburst are initially fitted with combined DBB and 
PL model components in XSPEC. Model fitted disk tempera¬ 
ture (T,„ in keV), power-law photon index (T), and flux con¬ 
tributions from two types of model components are obtained. 
We then refitted all the spectra with the current version of our 
TCAF model and from the fit, accretion flow parameters, such 
as disk rate (m,i), halo rate (/«/,), location of the shock (X s ) and 
compression ratio (R) are extracted. 

3.1. Spectral Data Fitted by TCAF Solution and by 
Combined DBB and PL model 

A combined conventional DBB and PL model fit in 2.5 - 
25 keV energy range RXTE PCA spectra provides us with a 
rough estimate of flux contributions originated from both ther¬ 
mal (from DBB) and non-thermal (from PL) processes around 
a BH. From this, we also get an idea about the evolution of the 
average temperature of the accretion disk and spectral states 
by monitoring variations of T,„ and T factors. However, from 
the variation of the TCAF fit parameters (such as two types ac¬ 
cretion rates, nij & nip, shock parameters, X s & R: and derived 
physical parameters, (such as the accretion rate ratio (ARR), 
shock temperature T s ^k, shock height li s hk, ratio between Iim 
to X s ), the accretion flow dynamics and geometry variation 
during the outburst phase become very evident. In Appendix 
Table I, all these fitted/derived parameters are written in a tab¬ 
ular form with estimated errors. 

Figures 1-3 show the variation of X-ray intensities, QPO 
frequencies along with TCAF and combined DBB, PL model 
fitted and derived (from TCAF) parameters. In Fig. la, vari¬ 
ation of the background subtracted RXTE PCA count rate in 
2-25 keV (0 - 58 channels) energy band with day (MJD) 
is shown. In Fig. lc, variation of TCAF fitted total accretion 
rates (combined Keplerian disk and sub-Keplerian halo rates) 
in the 2.5 - 25 keV energy band are shown. For comparison, 
combined DBB and PL model fitted total flux variation with 
day (MJD) is shown in Fig. lb. We observe that the variation 


of the TCAF fitted total flow rate (Fig. lc) is different from 
the flux variations in Figs, l(a-b) especially in early and late 
stages. In Fig. Id, variation of Accretion Rate Ratio (ARR, 
defined to be the ratio of sub-Keplerian halo rate ni/, and Kep¬ 
lerian disk rate mf) is plotted. Observed QPO frequencies (of 
only dominating primary QPOs) are shown in Fig. le. Dur¬ 
ing the entire phase of the current outburst, only three spec¬ 
tral classes, such as hard (HS), hard-intermediate (HIMS), 
and soft-intermediate (SIMS) are observed. Strangely, the 
soft state (SS) is not prominent and possibly missing. The 
sequence is found to be: HIMS (rising) —> SIMS —> HIMS 
(declining) —> HS. We believe that the absence of HS in the 
rising phase is due to observational constraints. The detailed 
behavior of these spectral states and other reports in the liter¬ 
ature on them are discussed in the next Section. 

In Fig. 2, variation of TCAF fitted and derived shock pa¬ 
rameters, together with combined DBB & PL model fitted re¬ 
sults are shown. In Figs. 2(a-b), variation of DBB temperature 
{Ti n in keV) and power-law photon index (T) with day (MJD) 
are shown. In Figs. 2(c-d), TCAF fitted shock location (X s 
in r g ) and compression ratio (R) are plotted with day. In Fig. 
2(e-f), variation of shock height ( h s hk in r g ) and temperature 
{T s hk in 10 10 K, which is the initial temperature of CT95 it¬ 
eration process), derived from X s & R and using Eqs. 4 & 5 
respectively of DMC15, are shown. In Fig. 2g, ratio between 
shock height and location are plotted. In Figs. 3(a-b), varia¬ 
tions of DBB flux from DBB & PL model fits and Keplerian 
disk rate from TCAF fits with day (MJD) are compared. Sim¬ 
ilarly, in Figs. 3(c-d), variations of PL flux and sub-Keplerian 
halo rate from these respective models are compared. Clearly 
there are ‘some’ similarities in each pair of compared quanti¬ 
ties, but not totally, since the power-law flux is a function of 
disk rate as well. In Figs. 4(a-c), TCAF fitted 2.5 - 25 keV 
background subtracted PCA spectra of three different spectral 
states (selected from approximate middle of each state to get 
better understanding, marked as a, b and c in Col. 1 of Ap¬ 
pendix Table I) along with residual^ 2 are shown. In Fig. 4d, 
we show unabsorbed theoretical spectra in 0.001 - 3950 keV 
energy range, which are used to fit observed spectra presented 
in Figs. 4(a-c). 

3.2. Evolution of Spectral and Temporal Properties during 
the Outburst 

Detailed temporal and spectral properties of this candi¬ 
date during this outburst are discussed by several authors on 
the basis of X-ray variability, QPO observations, spectral re¬ 
sults based on inbuilt XSPEC model fits, such as, power-law 
and multi-color DBB components, etc. (Munoz-Darias et al. 
2011; Kalamkar et al. 2011, etc.). However, since TCAF pro¬ 
vides us with variation of physical parameters from the spec¬ 
tral fits on a daily basis, it may be possible find a pattern to 
correlate with spectral classes. We found interesting correla¬ 
tion in 2010 outburst of H 1743-322 (MDC14) and in 2010-11 
outburst of GX 339-4 (DMC15). In Figs. Id and le, we see 
variations of ARR and QPO frequencies with time in MJD. 
We find that ARR, total flow/accretion rate (md+nih), shock 
locations, compression ratios, etc. in conjunction with QPOs 
provide a better understanding on the classification of spectral 
states. This will be discussed below. 

(i) Hard-Intermediate State in the Rising phase: RXTE has 
started observing the source three days after its discovery. 
Probably the initial low-hard state of the source in the ris¬ 
ing phase was missed. For the first 3 days of our obser- 
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vation (from MJD = 55467.19 to 55469.09), source was in 
a hard-intermediate state with increasing thermal DBB flux 
(also, Keplerian disk rate; see. Figs. 1 & 3). From Fig. 3d, 
it is seen that during this phase, TCAF fitted sub-Keplerian 
halo rate shows a rapid fall, although combined DBB & PL 
model fitted non-thermal PL flux is not changed significantly. 
This fact can be better understood by observing variations of 
ARR in Fig. Id. The shock moved rapidly (from ~ 354 to 
207 r g ) towards the black hole with reducing shock strength 
and height. This third observed day (2010 Sept. 30, MJD 
= 55469.09) is the transition day from hard-intermediate to 
soft-intermediate spectral states. QPO frequency increased 
monotonically (from 1.607 Hz to 2.723 Hz) with time (day) 
and ARR values decreased rapidly (from 3.259 to 0.590). 
According to propagating oscillatory shock (POS) model 
(Chakrabarti et al. 2005, 2008, 2009; Debnath et al. 2010, 
2013; used to explain mono tonic evolutions of QPO frequen¬ 
cies during rising and declining phases of the outburst), QPO 
rises rapidly till the compression ratio R reaches nearly around 
unity as the post-shock is cooled down. This is also seen in 
this outburst as well (Fig. 2d). 

(Hi) Soft-Intermediate State: 

The constancy of ARR lasted till the total rate as well as 
non-thermal (PL) flux or halo rate started rising suddenly on 
2010 Nov. 01 (MJD = 55501.23). This phase continued for 
~ 32 days, where sporadic QPOs are observed with very little 
changes in ARR, 7)„, T, R, T s hk , h s hk , and X s are observed. 
During this phase, the total X-ray intensity, flux or flow rate 
initially increased and then decreased mainly because of sim¬ 
ilar variations in thermal DBB flux or Keplerian disk rate. 
In this phase of the outburst, non-thermal PL flux shows de¬ 
creasing pattern, although sub-Keplerian halo rate initially de¬ 
creases, and then becomes more or less constant. On the soft- 
intermediate to declining hard-intermediate transition day, a 
rise in ARR value due to the effect of sudden rise in non- 
thermal PL flux/halo rate is observed. On this day, QPO fre¬ 
quency had a maximum (5.951 Hz). 

(vi) Hard-Intermediate State in the Declining phase: This 
state continued for the next ~ 3 days, starting from the SIMS- 
H1MS (declining) transition day. During this phase QPO fre¬ 
quency decreases rapidly from 5.951 Hz to 2.563 Hz. ARR is 
found to increase slowly with a rise in halo rate compared to 
the disk rate (see. Figs. 1 & 3). Rapid decrease in power-law 
photon-index (r) also indicates that spectrum start to become 
harder from the day one of this state. A slow movement of the 
receding shock with little increment in the compression ratios 
and shock heights are observed during this phase of the out¬ 
burst. Nov. 05, 2010 (MJD = 55504.06) is the transition day 
from declining hard-intermediate to hard state. Interestingly, 
on this day, ARR is locally maximum (=0.327) and QPO fre¬ 
quency starts to decrease slowly after that. Precisely this be¬ 
havior was seen in our earlier TCAF fits on other BHCs (see, 
MDC14, DCM15) as well. 

(vii) Hard State in the Declining phase: 

The source is observed in this spectral state till the end 
of the observation of 2010 outburst starting from the transi¬ 
tion day. In this state, ARR (from 0.327 to 0.278) as well as 
observed QPO frequencies (from 2.563 Hz to 1.638 Hz) de¬ 
crease monotonically as in other objects fitted by TCAF. A 
slow decrease in PCA count rate, total (DBB+PL) fluxes and 
total flow (disk+halo) rates are observed with a similar de¬ 
creasing trend in both thermal (DBB and mf) and non-thermal 
(PL or nif) flux/rate components. This is because the supply 
rate is dwindling after the peak outburst is over (see. Figs. 1 


& 3). A fast receding shock (from ~ 103 to 416 r g ) with rise 
in compression ratio and shock height are observed. At the 
same time, during this phase of the outburst, a decrease in 
power-law photon index is observed, which indicates that the 
spectrum becomes harder with a clear dominance by the sub- 
Keplerian halo and non-thermal power-law photons. Shock 
temperature ( T s /,k of initial iteration) values are found to de¬ 
crease monotonically with time (day). 

4. DISCUSSIONS AND CONCLUDING REMARKS 

We study the evolution of spectral properties of Galac¬ 
tic transient black hole candidate MAXI J1659-152 during 
its first (2010) X-ray outburst using the current version vO.l 
of two component advective flow (TCAF) solution based 
model after its inclusion as a local additive table model 
in HEASARC’s spectral analysis software package XSPEC 
(DCM14). This has been done with a model fits file using 
~ 4 x 10 5 theoretical spectra which are generated by vary¬ 
ing five model input parameters (two types of accretion i.e., 
Keplerian disk mj, sub-Keplerian halo mu rates; two types of 
shock parameters: location X s and compression ratio R: and 
the mass of the black hole Mbh) to the modified CT95 code 
(see, DMC15 for details). We re-fitted all these spectra of 
MAXI J1659-152 with combined DBB and PL model com¬ 
ponents to get a rough estimate about the variations of the 
thermal (DBB) and non-thermal (PL) fluxes during the out¬ 
burst and compare these with our TCAF fitted results (see. 
Figs. 1-3). In Appendix Table I, detailed results of our spec¬ 
tral fit with observed QPO frequencies are presented. 

The entire period of the 2010 outburst of MAXI J1659- 
152 appears to have three spectral classes: hard, hard- 
intermediate, and soft-intermediate. It did not reach the soft 
state. When we study variations of TCAF parameters in these 
states, we find that there is a pattern in how the rates, ARR, 
QPO frequency etc. behave. These behaviours are similar to 
what were reported in other sources, (see, MDC14, DMC15). 
Specifically, we see a local maximum of ARR during the tran¬ 
sition between hard-intermediate to hard states in all these 
sources. Clearly more objects need to be fitted before any 
firm conclusion can be drawn. 

It is interesting that unlike other sources (MDC14, 
DMC15), this object exhibited no soft-states during this 
outburst according to our model. Only for two days, 
MJD=55481.71 and MJD=55485.16 the fluxes are higher 
(Table I of Appendix), but observation of LFQPOs and the 
presence of a dip on MJD=55483.92 in between, suggests that 
the state is not soft, but soft-intermediate, van der Horst et al. 
(2013), using a disk irradiation model called DISKIR (with 
number of free parameters significantly higher than TCAF 
with irradiation from CENBOL on the Keplerian disk), sug¬ 
gests that the object might have gone to a soft state. However, 
the photon index (see. Fig. 5 of van der Horst et al. 2013) of 
those specific days showed significant error bars and thus it is 
uncertain if soft states were reached. It is also possible that 
the inclination angle might also have played a role in harden¬ 
ing the spectra (as discussed in another source by Motta et al. 
2010 and by explicit Monte-Carlo simulation by Ghosh et al., 
2011 ). 

Although low frequency (~ 0.01 - 30 Hz) QPOs are ob¬ 
served almost three decades ago, there is a debate on the ori¬ 
gin of this temporal behavior in Fourier transformed power- 
density spectrum of the X-ray intensity variation. Accord¬ 
ing to shock oscillation model (SOM) of Chakrabarti and his 
collaborators in mid-90s, it can occur when Compton cool- 



Accretion Dynamics of MAXI J1659-152 with TCAF 


5 


ing time scale roughly agrees with infall time scale (Molteni, 
Sponholtz & Chakrabarti, 1996) or due to non-satisfaction of 
Rankine-Hugoniot conditions to form a stable shock (Ryu, 
Chakrabarti & Molteni, 1997). Recent numerical simulations 
of Garain et al. (2014) also demonstrated this in presence of 
Comptonization. According to SOM, the QPO frequency is 
inversely proportional to the infall time (f,„/ a zz) in the post¬ 
shock region. It also has been observed that these QPO fre¬ 
quencies show monotonically increasing (during rising phase 
of the outburst) or decreasing (during declining phase of the 
outburst) nature in hard and hard-intermediate spectral states. 
This evolution of the QPO frequencies can be well fitted with 
the POS model, which is nothing but time varying form of the 
SOM. Movement of the shock inward could be due to rapid 
cooling and consequent collapse of the CENBOL (in the ris¬ 
ing phase) and outward, due to the lack of cooling in the de¬ 
clining phase (Mondal, Chakrabarti & Debnath, 2015). 

In soft-intermediate spectral states, sporadic QPOs are ob¬ 
served, which may be due to appearance or disappearance of 
the oscillating component, namely, CENBOL by intrusion of 
strong toroidal magnetic fields. Strong sporadic jets are also 
seen in these states (e.g., Nandi, et al. 2001; Radhika & Nandi 
2014). Our current TCAF model takes care of the combined 
effects of CENBOL and the outflow. Separation of the effects 
of CENBOL and jets is possible from more detailed modeling 
of timing properties and will be incorporated in a later version 
of TCAF 

Recently, it has been shown (DMC15, MDC14) using ex¬ 
amples of 2010-11 outburst of Galactic BHC GX 339-4 and 
2010 outburst of Galactic BHC H 1743-322 how spectral state 
transitions may be triggered when the relative ratio of the ac¬ 
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cretion rates, namely, ARR, vary in specific ways. The na¬ 
ture of variation of QPOs (when observed), shock locations, 
strengths etc. are also very specifically. Exactly same type 
of variation of the fitted parameters are also seen for the cur¬ 
rent source MAXI J1659-152. From the observed variation of 
Keplerian rates in these objects, we believe that an outburst 
is triggered due to a sudden rise in viscosity and is turned off 
due to the reduction in viscosity (CT95, Ebisawa, et al. 1996; 
Chakrabarti, Dutta & Pal, 2009). It is possible that this ob¬ 
ject belongs to a category, with short orbital period, where 
accretion disk around the black hole is mostly dominated by 
the wind accretion compared to disk accretion. The Keplerian 
disk is always immersed inside a strong sub-Keplerian halo. 
So, the soft state may be difficult to achieve. In future, we 
will make detailed spectral and temporal study of other such 
objects (for e.g., XTE J1118+480 of orbital period ~ 4.1 hrs, 
Gonzalez-Hernandez et al. 2013; Swift J1753.5-0127 of or¬ 
bital period ~ 3.2 hrs, Zurita et al. 2007) during their X- 
ray outbursts to check if flow dynamics of these sources also 
follow a similar trend. Prediction of QPO frequency from 
TCAF solution fitted shock parameters (X s & R; DCM14), 
and comparative study with POS model solution will be pub¬ 
lished elsewhere. 
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Fig. 1. — Variation of (a) 2 — 25 keV PC A count rates (cnts/sec), (b) combined disk black body (DBB) and power-law (PL) model fitted total spectral flux in 
2.5 — 25 keV range (in units of 10 -9 ergs cm~ 2 s _1 ), (c) TCAF model fitted total flow (accretion) rate (in M.Edd\ sum of Keplerian disk, rtid and sub-Keplerian 
halo riih rates) in the 2.5 — 25 keV energy band, and (d) Accretion Rate Ratio (ARR; ratio between halo and disk rates) with day (MJD) for the 2010 outburst of 
MAXI J1659-152 are shown. In the bottom panel (e), observed primary dominating QPO frequencies (in Hz) with day (MJD) are shown. The vertical dashed 
lines indicate the transitions of between different spectral states. 



Fig. 2.— Variation of combined DBB and PL model fitted (a) disk temperature Ti n (in keV), and (b) PL photon index (r) with day (MJD) are shown in top two 
panels. Variation of TCAF model fitted/derived shock (c) location (X s in r g ), (d) compression ratio ( R ), (e) temperature ( T s hk in 10 10 K), (f) height (h s hk in r g ), 
and (g) ratio between h s hk & X s , with day (MJD) are shown. The shock height, and temperature are derived from Eqs. 4 & 5 respectively of DMC15. 
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Fig. 3.— In top panel (a), the variation of combined disk black body (DBB) and power-law (PL) model fitted DBB spectral flux and in panel (c), the variation 
of PL spectral flux (both in units of 10 -9 ergs cm ~ 2 s -1 ) in 2.5 - 25 keV energy range are shown. In panel (b), the variation of TCAF model fitted Keplerian disk 
rate and in bottom panel (d), the variation of sub-Keplerian halo rate m/, (both in MEdd ) in the same energy band are shown. 
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Fig. 4.— TCAF model fitted 2.5 — 25 keV PCA spectral flux, in units of NormalizedCounts s 1 keV 1 with variation of Ax, selected from three different 
spectral states: soft-intermediate, hard-intermediate, and hard states respectively are shown in plots (a-c). The observation IDs (where Z=95118-01) and fitted 
values of x^. ed are written down in the plots. In bottom-right plot (d), unabsorbed TCAF model generated spectra, which are used to fit these spectra are shown. 
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Appendix I 

2.5-25 keV Combined DBB-PL and TCAF Model Fitted Parameters with QPOs 
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Here X=95358-01, Y=95108-01, and Z=95118-01 mean the initial part of the observation Ids, and (a-c) mark TCAF model fitted results for three different states, 
presented in Fig. 4. Intermediate void space mark state transitions from HIMS—>SIMS, SIMS—>HIMS, and HIMS—»HS respectively. 

Ti n , and T values indicate combined DBB and PL model fitted multi-color disk black body temperatures in keV and power-law photon indices respectively. 

* DBBf, PLf represent combined DBB and PL model fitted 2.5-25 keV fluxes for DBB and PL model components respectively in units of 10 -9 ergs cm~ 2 s~ l . 
nth, and raj represent TCAF fitted sub-Keplerian (halo) and Keplerian (disk) rates in Eddington rate respectively. X s (in Schwarzchild radius r g ), and 
R are the model fitted shock location and compression ratio values respectively. h s hk (in r g ) and T s hk (in 10 10 K) are the shock height and temperature values 
derived from Eqs. 4 & 5 of DMC15 respectively. 1 ’ Here, frequencies of the principal QPO in Hz are presented. DOF means degrees of freedom of the model fit. 
Note: average values of 90% confidence ± error values obtained using ‘err’ task in XSPEC, are mention as superscript of the spectral fitted/derived parameters. 








